Electronic properties of the pseudogap system 
(TaSe4)2l 

Nic Shannon and Robert Joynt 

Department of Physics, University of Wisconsin-Madison, 1150 University Avenue, 
Madison, WI 53706, USA 

Abstract. 

The room temperature "metalhc" properties of the quasi-one-dimensional material 
(TaSe4)2l differ markedly from those expected of either a Fermi or a Luttinger Liquid, 
showing strong signs of a suppression of the density of states at the Fermi level. We 
present evidence for the existence of strong quasi-static fluctuations of structural 
order with long correlation length. These fluctuations produce a pseudogap in the 
density of states. We compute the temperature dependence of the optical and DC 
conductivities of (TaSe4)2l in its conducting phase, the nature of its core hole spectra, 
the NMR Knight Shift and relaxation rate. Predictions for these quantities are made 
on the basis of a Lee, Rice and Anderson model. This model represents the simplest 
theory of a pseudogap, and gives satisfactory agreement with experiment in the cases 
where comparisons can be made. In contrast, the predictions of a strongly correlated 
(Luttinger Liquid) model appear to to contradict the data. The chief remaining 
discrepancy is that the gap deduced from transport quantities is less than that observed 
in photoemission. We discuss some possibilities for resolving this issue. 
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1. Introduction 

(TaSe4)2l is among the most widely studied of quasi one-dimensional materials. First 
synthesized in 1984 it provides a good example of a "metallic" system undergoing 
a Peierls transition to a nearly commensurate charge density wave ground state. It is 
unusual among low dimensional systems in the uncommonly high temperature (263 K) 
at which it becomes ordered, and in the way in which that ordering affects its room 
temperature properties. It has attracted interest as a Luttinger liquid candidate, based 
on its photoemission spectrum We argued in an earlier paper [|] that none of 

the observations on (TaSe4)2l are really consistent with the Luttinger picture, and that 
the photoemission is best explained by its proximity to a charge density wave state — 
the energy scale associated with the gap which opens in the charge density wave state is 
clearly visible in the room temperature properties of (TaSe4)2l, as was shown by early 
measurements of magnetic susceptibility p. (TaSei)2l is clearly a 'pseudogap' system. 

Below 263 K, the resistivity of (TaSe4)2l shows an activated behavior, and 
measurements of its optical conductivity and of the dispersing features measured in 
angle-resolved photoemission (ARPES) spectra also give evidence of a gap. Strangely 
however, the gap (or, above 263 K the pseudogap) measured by different experiments 
is not of the same size, as is illustrated in Table 1, where gaps and their corresponding 
mean-field transition transition temperatures are listed. 



Table 1. Estimates of gap size and transition temperatures. 





A 


T3D 


7m F 


DC Conductivity (e.g. @) 


180 meV/2100 K 


263 K 




Optical Conductivity |^ 


200 meV/2300 K 






ARPES § 1 


520 meV/6000 K 




892 K 


Magnetic Susceptibility Q 






860 K 



Roughly speaking, the gap in ARPES spectra is 2 to 3 times that in other 
experiments. We discuss below how this discrepancy may be resolved. 

Coincidentally, observation of the pseudogap state in underdoped high-Tc 
superconductors has caused great excitement. The nature of this state is presently 
intensely debated, with no consensus having been reached. In this context, it is 
important to look in some detail at (TaSe4)2l which, we believe, is a much simpler 
pseudogap material. In fact we will contend that its properties, with the single exception 
of the gap discrepancy, can be understood in a straightforward, and even rather crude, 
theory. 

The dispersing features measured by ARPES in the conducting phase were found 
to be quite well described by a simple Lee, Rice and Anderson (LRA) model in which 
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the fluctuations of lattice order associated with the Peierls transition are viewed as a 
temperature dependant ensemble of static potentials in which the electrons move 0, 
and electron-electron interaction is explicitly not included. 

In this paper we extend the analysis begun in |^ to optical and DC conductivities, 
making predictions for the temperature dependence of these quantities in the conducting 
phase on the basis of the same simple LRA model. We will also discuss what might be 
learned from NMR and core hole photoemission experiments on (TaSe4)2l, considering in 
particular what might be established about spin-charge separation from the comparison 
of the two. In all cases we compare these predictions both with existing experimental 
data and with the behavior expected of models based on strong electron-electron 
interaction. For greater readability, ancillary technical details are given in appendices. 

2. The Lee, Rice and Anderson model 

Conduction in (TaSe4)2l takes place along one-dimensional chains of Ta atoms |l|. The 
material undergoes a Peierls transition at = 263K to a state where the lattice is 
distorted by a condensed transverse acoustic phonon mode. This may be thought of 
in simple terms as a slightly incommensurate tetramerization of the tantalum atoms 
perpendicular to the chain axis. Superlattice reflections appear in X-Ray spectra at 
q = (±0.05, ±0.05, ±1.085), confirming incommensurate charge density wave order not 
quite aligned with the chain axis [0]. In the chain, this corresponds very nearly to a 2kF 
fluctuation. 

The proper model for what is then a one-dimensional electron-phonon problem is 
given by the Frohlich Hamiltonian 



where 




and cl and 6j are (respectively) creation operators for electrons and phonons with 
dispersion e{k) and uj{q). Uq is the Fourier transform of the lattice displacement, and 
g{q) the electron-ion coupling. Since electron spin enters into the problem only in 
appropriate factors of two, it will be suppressed in our notation. 

The physics of this Hamiltonian has been widely studied for many years and is 
highly nontrivial. It is well known both that the one dimensional lattice is unstable 
against distortion, and that thermodynamic fluctuations prevent a transition to a state 
with long range order from occurring at finite temperature in any truly one-dimensional 
system. 
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The possibility of phonons with wave number Q = 2kp decaying into zero energy 
particle-hole pairs leads to softening of the phonon spectrum for Q ^ 2kF, and at a mean 
field level the system described by the one-dimensional Frohlich Hamiltonian Eq. (1) 
undergoes a transition to a fully ordered charge density wave state at the temperature 
Tmf for which the frequency of the Q = 2kF phonon mode goes to zero p. This 
temperature is determined by the strength of electron-phonon coupling g{2kF), and is 
of the same order as the gap developed in the charge density wave state at absolute 
zero. Within the set of approximations usual for BCS theory one has 2A = 3.5/cb^mf 
but for real materials this relation is seldom exact. The three-dimensional charge density 
wave transition which occurs in (TaSe4)2l and other similar compounds is stabilized by 
interaction between different metallic chains, but occurs at a temperature much smaller 
than the gap energy, as can be seen in Table 1 where Tmf is certainly well above the 
exerimental transition temperatue. 

In the spirit of LRA, we identified the transition temperature = T^u = 263K 
with a crossover from a three dimensional Peierls-distorted mean field charge density 
wave state described by 

Hmf = E <k)clck + E[^*4-2fc,c, + Acl,^,,^c_k] (2) 

k k 

g{2kp){u2kp)- (3) 



to a state in which there are essentially uncorrelated fluctuations of charge density wave 
order on individual chains. We make no attempt to accurately describe the way in 
which this dimensional crossover takes place. Fluctuation effects above T^d are taken 
into account in the simplest possible way consistent with a mean-field ground state, 
replacing the order parameter for the Peierls distortion of the lattice A = -^g{2kp){u2kp) 
by a static external field which has a non-vanishing expectation value for other momenta 
not exactly equal to 2kF: 

HLRA=J2<f'Hck+ E [^q4'-q% + ^qc1,,+qC_,,] (4) 

k Q,k'>0 

vI/Q =-L^(g)(«(g)). (5) 

The classical field belongs to a thermal ensemble of potentials characterized by a 
mean sqaure "gap" scale {ip'^iT)) reflecting the size of fluctuations of lattice disorder, 
and an inverse coherence length ^~^(T). 

Since the main purpose of this paper is to explore the LRA model as a 
phenomenology for (TaSe4)2l we consign the derivation of the single-particle Green's 
function for electrons moving in this static phonon field to [Appendix A| , and reproduce 
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here only the result 

where the choice of sign ±ivfC,~^{T) is made according to whether we continue to the 
upper or lower half plane. For convenience we have set h = 1. 

Determining the temperature dependence of {ip'^{T)) and ^~^(T) forms an essentially 
independent problem; for a fully self-consistent phemonemology they should be found 
from experiment. The scale of (V'^(T)) is set by the size of the mean gap at T = 0, and 
it varies little over the range of temperatures in which we can apply the LRA theory 
to experiments on (TaSe4)2l, i. e. 263 < T < A30K. In contrast ^^^(T) varies quite 
strongly with temperature, and becomes extremely small at the transition temperature 



T3D. For all experimentally accesable temperatures VfC,~^{T) << {ip'^{T)). Our 
parameterization of the model is discussed in more detail in [Appendix B . 
In what follows we will use a spectral representation of the electrons 

All the anomalous features of the room temperature ARPES data are present in this 
spectral function. In particular it predicts the broad dispersing features, and "quasigap" 
structure observed by Terassi et al 0. 

It is important to note that the results of this paper assume the validity of (7) and 
nothing further, since all experimental quantities can be expressed as integrals over the 
spectral functions. Since this equation can be directly compared with experiment, we 
may regard our theory in two different ways. On the one hand, it may be viewed as a 
phenomenology in which the spectral function is taken from one experiment (ARPES) 
and used to develop a picture of several experiments, or it may be viewed as a direct 
test of the LRA model against experiment, without prejudice as to the suitability of 
LRA as a solution to the Frohlich Hamiltonian. 



3. Previous Work 

For completeness, we review previous work on the application of these equations to 
explain experimental data. An example of a fit to a room temperature ARPES spectrum 
for (TaSe4)2l made on this basis in Ref. |Q is shown in Figure]^. A fit to the temperature 
dependence of the uniform magnetic susceptibility of (TaSe4)2l by Johnston et al. 
is given in Figure ^. In both cases, the agreement of theory and experiment is excellent. 
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4. DC Resistivity 

The DC resistivity p(T) of (TaSe4)2l is usually presented in an Arrhenius plot so as 
to extract an ordering temperature Tjd and an activation energy Aq for the charge 
density wave state. For temperatures above T^d the data plotted in this form show no 
upturn — in fact the resistivity becomes nearly temperature independent, decreasing 



slowly over all higher measured temperatures |110|. This contradicts our expectations 
of an ordinary metal, where most scattering mechanisms increase in effectiveness with 
increasing temperature, and an upturn in p(T) is expected on the closing of the gap. 

In the LRA liquid we anticipate that scattering of electrons from fluctuations of 
charge density order will have two effects on conductivity: the suppression of the 
density of states at the Fermi energy, and the imposition of a finite lifetime on electrons 
propagating in momentum eigenstates along the chains. The former is reflected in 
the pseudogap visible in photoemission spectra and the spectral function Equation (7), 
the latter in the large temperature-dependent width of the dispersing features in the 
spectral functions. The temperature dependence of the conductivity of will therefore 
depend on the interplay between the gap and lifetime effects as parametrized by (ip'^lT)) 
and rH^)- 

We calculate the intrinsic DC conductivity of an LRA liquid directly from a current- 
current correlation function. The results of the usual Kubo formalism, evaluated to first 
order, may be expressed directly in terms of the spectral function of the system under 
consideration as : 



k 



We evaluate the sum over as a contour integral and perform the final integral over 
uj numerically. We do not include disorder in our calculations in this paper, and 
the DC conductivity is of course sensitive to disorder. We therefore add a small 
temperature independent part to the inverse coherence length ^~^(T), which we chose 
to be ^~^{300K). This affects the temperature dependence of the conductivity only 
near T3D. In order to present consistent results we will include this correction in later 
calculations of optical conductivities. 

Theory and experiment are shown in the interesting range from 26AK to 430K in Fig. 
^. There is qualitative agreement, in that the unusual slow decrease with temperature 
is reproduced. However, it is clear that a proper treatment of the crossover to three 
dimensions is needed to fit the data near T^d- 

To understand the qualitatively behavior of p(T), we can make an estimate of the 
scale of the temperature dependence somewhat in the spirit of the Drude model : 

a(T) ~ Z^LRA(e/) X t^^a, (9) 
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where Tlra is the effective hfetime of the electrons and Z^lraIc/), the density of states 
at the Fermi energy, has been chosen as a representative measure of the degree to which 
the gap has filled. 

An expression for Di^^a{uj) can be found analytically from (7); for = e/ it has the 
simple form 

Naively, we might expect 1 / Tlra to be given by the imaginary part of the self energy of 
an electron at k = kf, u = 0, but this diverges for ^ ^ 0, while the reciprocal of the real 



electron lifetime should tend to zero. However in the limit where f/^~^(T) / y {ip'^iT)) 
it is possible to show that the spectral function of Eq. (7)) can be rewritten as the 
sum of two Lorentzians, with Vf^~^{T)/2 playing the role of a lifetime (peak width). A 
naive prediction for the DC conductivity would then be : 

Since {tp'^{T)) decreases only slightly over experimentally relevant temperatures, this 
represents a weakly temperature-dependent increasing conductivity, as is observed. But 
note that the conductivity we calculate is anything but temperature-independent over 
the remainder of the (greater) range of temperatures which we may access in our model. 
By fitting the numerical results we find it is rather well described by a high power law 
in the reduced temperature (T/TufY, x ~ 10, with the weak temperature dependence 
of the conductivity at low temperatures being due to its suppression by many factors of 
T/Tmf- 

Thus the very weakly temperature dependent DC conductivity is well explained 
by a cancellation of density-of-states and lifetime effects. This fails to capture the 
interesting step-like feature seen at the transition temperature in an Arrhenius plot. 
This is not surprising, as the dimensional crossover is not included in a realistic fashion 
in our theory. This (near) independence of DC conductivity on temperature over so 
large a range of temperatures is sufficiently unusual that we consider it to be a striking 
confirmation of the LRA theory. 



5. Optical conductivity of the LRA Liquid 

The optical conductivity of (TaSe4)2l has been measured over a range of temperatures 
IP and also shows a surprising lack of evolution with temperature. 

A truly gapped state, such as that described by Eq. (2) cannot support charge 
carrying excitations (or indeed any excitations) at energies less than the 2A needed 
to promote charge across the gap. The measured optical conductivity should therefore 
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"turn on" abruptly at ^ 2A. Since the density of states diverges immediately above 
and below the gap, we expect a very pronounced asymmetric peak at u = 2A with 
almost all weight on the high energy side. 

Such a peak is indeed observed in the charge density wave state of (TaSe4)2l. 
However, it is rather less asymmetric than might have been expected. Away from the 
maximum the spectra grow as a;^ on the low- and fall away as l/u^ on the high-energy 
side. Since we are only interested in the excitations from the band of electrons which 
carry currents in (TaSe4)2l, we will ignore the presumably interband contribution which 
sets in beyond 1.2eV. 

In a mean-field-like picture, we would expect the feature at uj ^ 2A ^ AOOmeV to 
become weaker, broader, and migrate to lower energies as the temperature is raised and 
the gap closes. The first two of these expectations are fulfilled, but over the range of 
temperatures 15K — 4007^^ the position of the maximum changes by a few meV only, 
and not by the ~ 120meV that might have been expected from our parameterization of 
the model. 

The only qualitative change in the peak as (TaSe4)2l is heated through T3D is the 
loss of the simple uo"^ behavior on the low energy side — on a log-log scale the peak 
grows a small shoulder at a litUe below A. On the high energy side the peak still falls 

away as 

In fact the Hamiltonian in Eq. 2 is not adequate by itself to describe the dynamics 
of the ordered state. The lattice distortion may vary slowly in space and time p3[] , 
leading to low energy collective modes in the charge density wave which are observed at 
finite frequency in Blue Bronze even in its conducting phase . Since we are mostly 
concerned with pseudogap effects in the single-electron properties of (TaSe4)2l however, 
we will not consider such collective excitations here. 

Under this approximation, the optical conductivity of the Lee, Rice and Anderson 
model can be calculated directly from a Kubo formula expressed in terms of electron 
Green's functions, in the same way as the DC conductivity. In this case we start from 

e'^ de 



2m? J-00 2tx 



nf[e + io) — Uj 



^ J k 



Y,k^Mk,e)A{k,e + uj) (12) 



which obviously has Eq. (8) has its a; — > limit. Once again we perform the sum over 
k analytically and the integration over u numerically. Typical results are presented on 
a log-log scale in Figure ^. They may be characterized by 

• Lorentzian peak at a; = with width of the order the naive quasiparticle inverse- 
lifetime vf^"\T). 

• Strongly suppressed response over the range of energies u ^ vfC,^^{T) — >• 



Shoulder at uj ^ J {ip'^iT)) — v'p^^^iTY (for temperatures sufficiently close to Tgo). 
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• Sharp asymmetric peak at ~ y (ip'^iT)) — Vp^~^(Ty with temperature dependent 
position, width and height. 

• Asymptotic decay as l/u^ at high frequencies. 

These features of the spectra and their temperature dependences are again in good 
quahtative agreement with those taken on (TaSe4)2l, although it is important to note 
that we have been forced to abandon the value of the zero temperature gap Aq = 0.52eV 
required to fit photoemission in favour of the smaller value Aq = 0.2eV suggested by 
transport measurements. This value is listed in Table 1. We have explicitly checked 
that the — > limit of the predicted optical conducivity matches our prediction for the 
DC conductivity, providing a useful self-consistency check on each calculation. 

The important differences been model and experimental spectra are 

• Greater temperature dependence of the peak position. 

• Lack of weight within the gap — experimental spectra have more weight within 
the "gap" especially at low temperatures. 

Interestingly recent higher resolution photoemission spectra for (TaSe4)2l also show 
less temperature dependance in the peak position than the earlier data which we fitted 



in terms of the LRA model |]T5|. Without attempting a detailed analysis, we note 
that optical conductivity spectra taken on Blue Bronze are qualitatively very similar 
to those taken on (TaSe4)2l. Spectra for the family of organic quasi one- dimensional 



systems {TMTSF)2X, on the other hand, show important differences in detail JIB 



We wish to stress that measurement of the different gap scales by different 
experiments, and the excess of weight in the gap are problems not of the conducting 
phase only, but also need to be addressed in the ordered phase. To do this it may be 
necessary a) to go beyond the approximation of static lattice distortion which underlies 
both the mean field theory of the ordered state and the LRA model and b) to include 
the effects of electron-electron interaction and/or impurities in a more careful way. 

One way of obtaining a better fit to optical conductivty spectra in the ordered phase 
would be to convolute the expected mean field response function (which is singular at 
2 A) with a Lorentzian "damping function" Jl^. This however leads to the transfer of 
rather too much weight into the subgap region [uj < 2A). It has been argued that the 
effect of both thermal and quantum fluctuations in the ordered phase can be modeled 
with a zero mean "white noise" potential (i.e. one with only delta function correlations in 
real space) |jl8[ . Predictions based on this model have been sucessfully applied to optical 
conductivity spectra for the charge density wave states in KCP {K2Pt{CN)^BrQ ^?)H20) 
|T9| and Blue Bronze {KMnO^) [0. In the case of (TaSe4)2l 0], while the spectra 



showed some deviations from the predicted behavior, the temperature dependance of 
spectra for T < Tg^, oj < 2 A led the authors to conculde thermal that lattice fluctautions 
are indeed an important effect, even in the fully ordered charge density wave state. 
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6. Comparison with Strongly Correlated Models 



The properties of (quasi) one-dimensional strongly correlated electron systems have 
been studied theoretically for more than twenty years and good reviews exist [^. We 
shall only highlight the transport properties of interacting systems which differ radically 
from non-interacting ones in ways which are relevant for the experiments discussed 
above. 

We may characterize the Luttinger Liquid (LL) state as a) spin-charge separated - 
low energy spin and charge excitations are independent of one another; b) gapless — 
low spin and charge excitations are bosonic collective modes with "acoustic" dispersion 
^p{a) = Vp[a)\(l\', and c) critical - all correlation functions decay with power laws 
determined by two parameters interaction Kp and K„. (We assume spin rotation 
symmetry below, in which case K^r = 1. For a recent review of LL physics, see p2|). 

These properties are most conveniently illustrated by the asymptotic form of the 
single particle correlation function 



G{x,t) 



1 

2^ 



Aktx 



+ 



g tic j-' X 



— Vpt^Jx — Vfjt ^X + Vpt^Jx + Vfjt 



x^ 



where x,Vpt,v^t ^ A, a cutoff (interaction range), and 



e 



-2 



(13) 



(14) 



Local repulsive interactions lead to 1/2 < i^p < 1 [0 < < l/8v^]. 

It is immediately evident from this correlation function that the Luttinger Liquid 
does not support Fermionic quasiparticles — it has no pole structure in the propagator 
and therefore no Fermi surface. That the density of states vanishes as a power law 



Diuj) 



.20 



(15) 



and that the Fermi surface "step" is replaced by a weaker algebraic decay \n{k) — 1/2| ~ 
\k — kfl"^^ follows from setting x = (t = 0) in (13) and Fourier transforming on t{x). In 
fact the spectral properties of the models leading to (13) have been worked out in some 
detail [^, and represent a convolution of the responses of the spin and charge sectors. 
For a general "universal" parametrization of the interaction, the single electron spectral 
function has two singularities, dispersing as Vp{k — kj) and Va{k — kf). The extent 
to which the spin-charge separated model is not universal in its spectral properties is 



dicussed in 24 



Like the free electron gas, the Luttinger liquid is a perfect conductor, and the finite 
conductivity of any experimental realization of a Luttinger liquid would be determined 
by the effect of additional terms in the Hamiltonian arising from disorder or coupling 
to external fields (crystal field, phonons, etc.). These are strongly renormalized by 
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interaction and scaling arguments lead to power law temperature dependence of 
(Jo, and power law behavior in uj for (j{uj) arising from those (Umklapp) terms which 
open a gap to charge excitations ||21|, |26[] . 

p(T) ~ T'-" (16) 
a{uj) ~ uj-'' (17) 
z/ = An^Kp - 5 (18) 

where n is the commensurability of the system (inverse of number of electrons per site). 
(TaSe4)2l is naively a quarter filled system, in which case n = 2 and 

u = 16Kp - 5 (19) 

Room temperature optical conductivity and ARPES measurements show clear evidence 
of a gap scale A, which means that the conducting phase of (TaSe4)2l cannot be a LL 
by our definitions above. Nonetheless we might can still look for evidence of strong 
correlation effects for u 3> 2 A, with the understanding that structure at lower energies 
will be non-universal and rather more complicated. 

Fitting the uj~^ tail of the optical conductivity leads to the value Kp = 1/2 - the 
extreme limit of what may be accomplished with local electron-electron interaction, but 
not incompatible with strong electron-phonon interaction. In this case we expect a gap 
to open in the spin sector and a renormalization 



= J^Kp (20) 
V m* 

where we have followed the convention usual in CDW literature of associating 
a renormalized electron mass m* with the dynamics of the charge density wave. 
Empirically this may be as much as a few hundred times m. 

Self consistently using the value Kp = 1/2 in (16) we would then anticipate 
(Jo ~ 1/T, which is entirely incompatible with the monotonically increasing conductivity 
measured experimentally. To attempt to extract a value of Kp by fitting the vanishing 
of the density of states at the chemical potential with the form 15 is clearly perverse in 
a system with so large a gap, and leads to even more extreme values of Kp ^1/3. 

It is worthwhile noting that the conductivity calculated above is the the 
hydrodynamic conductivity of the hquid (the response of an isolated system in 
equilibrium to an applied field), and not an experimental conductance. Predictions 
for the Landauer (two terminal) conductance of an isolated LL depend on how one 
considers it to be coupled to the current carrying leads of the experimental apparatus 
27| , p8| . As far as we are aware, the question of how a real system of many chains couples 



to external current carrying leads has not been considered in detail by any author. 

The LL is thus clearly a poor candidate for describing the low energy properties of 
the conducting phase of (TaSe4)2l: experiment shows clear evidence of a pseudogap at 
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room temperature in both spin and charge channels, and even above this gap energy 
scale it does not obey the universal properties scaling expected of a such a model. 

Another paradigm for a one-dimensional conducting state is the Luther-Emery (LE) 
Liquid, which is spin-charge separated and has the same charge excitations as a LL, but 
a gap to spin excitations. Moreover, while a LL supports both charge density wave and 
spin density wave excitations decaying with the same power law, the spin gap in the LE 
Liquid singles out charge density wave fluctuations, making it a natural candidate for 
describing quasi one-dimensional charge density wave systems. 

The LE Liquid should respond in the same way as the LL in all responses which 
couple to charge alone, and so it will have the same ideal hydrodynamic conductivity 
as a LL or a free electron gas. An important consequence of this is that its optical 
conductivity should obey the same power law behavior as that of a LL. The spin-gap is 
felt however, whenever a probe measures spin excitations (susceptibility measurements) 
and also if it couples directly to electrons, as in photoemission, tunneling experiments, 
and two terminal conductances, so all of these should show evidence of a (quasi-)gap. 
As both charge and spin sensitive probes reveal clear evidence of a gap in (TaSe4)2l, we 
do not feel that there are any strong arguments in favour of the LE Liquid in this case. 
That the LE Liquid may provide a rather better description of Blue Bronze has been 
suggested by Voit . 



7. Core Hole Spectra 

Core level X-Ray photoemission spectroscopy (XPS) provides a useful probe of the 
structure of a material, and of the low energy excitations of any charge carriers present 
in it. Generically these spectra have the form of a set of lines in an incoherent 
background, at energies many electron volts below the chemical potential. Because 
the wavefunctions of the electronic core levels are physically very small, they provide 
an essentially local probe of the structure of the material, with shifts in individual lines 
providing information about changes in the electric field and charge susceptibility at the 
sites of individual atoms. 

The details of the lineshape of electrons emitted from core levels also depend 
sensitively on the low energy excitations of any free charge. These excitations are 
in turn dominated by many-body effects and the asymmetries observed in lineshapes 
for metals are closely related to another much studied problem - the "Fermi Edge 
Singularity" (FES) in X-ray absorption. We will argue that, taken in conjunction with 
a complementary local probe of low energy spin excitations such as NMR, XPS (or, with 
some modification to our arguments, the FES) offers a novel means of distinguishing 
between different metallic and conducting phases. 

The calculation of many body effects in core level and X-ray edge response in three 
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dimensions for a Fermi Liquid has been extensively studied in a body of work commonly 
referred to as MND (Mahan-Nozieres-Dominicis) theory . 

In the case of ordinary metals it is enough to consider the noninteracting problem. 

H =Ho + V (21) 
Ho= J2<k)clck + e,,dU (22) 

k 

y = iT.yA+,c,dd^ (23) 

k,q 

where c|, creates a conduction electron with energy e{k) and creates a core electron 
with binding energy eh- In higher dimensions {k, q} are vectors {k, q}. Spin indices have 
been suppressed for compactness. We make the usual assumption that lineshape for the 
core level measured in an XPS experiment is proportional to the Fourier Transform of 
the (retarded) core hole Green's function 

Ghit)= -t9it){\d^it)d\) (24) 

broadened by an appropriate factor to allow for experimental resolution. In practise to 
extract a usable lineshape allowance must also be made for the finite lifetime of the core 
hole (modeled as a Lorentzian with width set by Auger decay processes) and in some 
cases a symmetric (gaussian) broadening of the line due to interaction with phonons. 
These issues are discussed in the associated literature USD . 



An important many body effect comes into play in the evaluation of this correlation 
function. As originally observed by Anderson [^], the ground states of a free electron 
system with and without the impurity are orthogonal, their overlap vanishing as 1/A^"°, 
where 

"o = iEf^V (25) 



is the total number of particles and Si is the phase shift in the l-th scattering channel. 
In three dimensions / has the natural interpretation of an angular momentum quantum 
number for scattering states. 

In metallic systems this orthogonality makes itself felt in XPS line asymmetries and 
in a contribution to the FES in X-Ray absorption; the number of zero energy particle 
hole pairs which may be made at the Fermi edge is essentially only bounded by the 
level spacing of the system, which in the thermodynamic limit, for a (i-dimensional 
system in volume L'^ vanishes as 1/L. An important consequence of this for XPS is 
that in removing a core electron from a metal, a very large number of conduction band 
electron-hole pairs are excited by the unscreened core hole. This converts the lineshape 
measured by and core hole photoemission experiment from a simple 6 functional form 

Gh{uj) oc 5{uj - cut) (26) 
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to a power law asymmetry 



Ciu) (X , . (27) 



T) 



where ujt and uj'rp are the threshold energies for emission of a core electron (naively given 
by eh). 

For free electrons the XPS exponent a differs simply by a factor of two from the 
orthogonality exponent for the overlap defined by Anderson. Within the Born 
approximation, 5i will be proportional to the modulus of the matrix element for the 
core hole potential in that channel: 5i ~ In one dimension there are only two 

scattering channels "forward" (g ~ 0) and "backward" {q ~ 2A;/), so, 

«'^^ = ^o(|VoP + |V^2./P) (28) 

where Uq = l/2nvf is the density of states at the Fermi surface. 

X-Ray spectroscopy of this type is a probe of the low energy charge excitations at 
a given lattice site, and as such essentially measures its local dielectric constant. We 
can infact generalize the result Eq. (28) to an arbitrary many electron system, using the 
relation 

a= limElV^.P^^^^^ (29) 

which we derive in Appendix C, provided that certain conditions on Xpil^^) K 
hold. 

There is also a closely related shift in the XPS line position 

AE= ElKr3^{Xp(g,0)} (30) 

Both of these parameters acquire temperature dependence in a system where the charge 
susceptibility changes with temperature. 

Returning to the case in point — (TaSe4)2l — since the density of states at the Fermi 
energy is zero in a mean field charge density wave state, and strongly suppressed for 
an LRA Liquid Eq. (10), we clearly expect very little asymmetry in XPS spectra for 
(TaSe4)2l in either its metallic or conducting phase. Such asymmetry as exists would 
follow the form in Figure |^. 

The temperature evolution of AE{T) may well be measurable. However since 
processes at all energy scales contribute to the real part of the local susceptibility, and 
the absolute value of the shift is cutoff-dependent there is no simple universal formula 
for the shift valid at all teperatures. The temperature dependance of a relative line shift 
must be found numerically. Results for this quantity are plotted for our parameterization 
of the model in Figure |^ 

This lack of asymmetry is to be expected: many body effects related to the existence 
of a Fermi surface will not occur in a system with a gap (or quasigap). However, such 
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an argument can only be applied in a system where the quasi-particles are in some 
sense electron — like. It is clearly unreliable in the case of a Luttinger Liquid where the 
quasi-particles are stable but are not electron-like in character. We consider this case 
below. 

The X-ray response of strongly correlated one dimensional metal (a Luttinger 
Liquid) has become something of a cause celebre in recent years. "Forward scattering" 
from a core level may easily be treated within the bosonization scheme used to 
diagonalize the interacting problem, leading to a suppression (enhancement) of the 
asymmetry measured in a system with repulsive (attractive) interaction between 



electrons, compared with that measured for non-interacting electrons |32 



In [Appendix E| we present an exact solution of the forward scattering contribution 
to XPS in a Luttinger liquid; here we simply reproduce the result in convenient form 

afW) = K.ia^ (31) 

P 

where Vp and Kp are the parameters of the charge sector of the LL, defined in 
Equation (13), and a^^ is the result for non interacting electrons within the Born 
Approximat ion . 

Unfortunately "backward scattering" introduces terms nonlinear in bose fields 
into the Hamiltonian and make the problem very much harder to treat analytically. 
Renormalization group analysis in fact suggests that strength of this nonlinear term 
flows to infinity, "breaking" the one-dimensional chain of atoms . Within this "open 
chain" interpretation of the RG results, the contribution to the orthogonality exponent 
ao from backward scattering is expected to take on the universal value of 1/16, and the 
orthogonality exponent has the value 

This result may be derived very elegantly from the use of boundary conformal field 
theory |Q . 



Equivalent results for XPS response have been derived in expansions about the 
"open chain" fixed point Generically these predict an exponent 

aE = '^Kp^a^^ + \ (33) 
vj •' 8 

at threshold, with a crossover to a different non-universal scaling law at a finite energy 
determined by the strength of the core-hole conduction electron coupling. A particularly 
elegant approach to the XPS problem which reproduces these and many other standard 
results is provided by the use of boundary conformal field theory |3^ . 

It should also be noted that not all authors have found their results in agreement 
with the "open chain" interpretation set out above, with one asserting that backward 
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scattering terms can lead to an enhancement of the FES A recent numerical 

treatment of the problem with references to earlier analytical work is provided by |36|. 

No predictions have yet been made for the temperature dependence of the XPS 
exponent in a LL, but we can speculate that these will also obey non-integer power 
laws controlled by Kp. As stated above XPS measures the leading u dependence of 
the local charge susceptibility a system. For a LL, we may divide the susceptibility 
into a regular part at Xpil ~ 0) and an anomalous part at Xp(g ~ 2kf). The former 
has the same structure as in a free electron gas (up to charge velocities); from scaling 
arguments we anticipate that the temperature dependence of the latter will be a power 
law controlled by the interaction parameter Kp. 

As stressed above, XPS is a zero energy probe and the scaling of Xp will be cut off 
by any gap. 

Core level studies for materials with a charge density wave groundstate (see, e. g. 
, Ref. |3^) usually reveal gross differences between spectra taken above and below the 
transition temperature. There are two fundamental reasons for these differences - the 
change in the local environment of core states due to a structural transition or charge 
density wave, and changes in the charge suspcetibilty induceded by the charge density 
wave transition, which modify the lineshape and position, as discussed above. 

In a one-dimensional tetramerized charge density wave state (assuming perfect 
commensurability) there are two inequivalent sites in each unit cell. The unbalanced 
charge accumulated due to the charge density wave will be different at each of these 
sites leading to a splitting in the XPS lines proportional to the magnitude of the charge 
density wave {p{2kf,T)), and therefore following it in temperature dependence. It is 
unlikely that such a splitting would be observed in an incommensurate charge density 
wave system hke (TaSe4)2l although this has been a subject of some contention ||3^ . 

Incommensurability can be expected to lead to the smearing of the core level 
threshold energy over a range of energies again proportional to the magnitude of the 
charge density wave, and therefore to lines whose width follows {p{2kf,T)). In the 
quasi-static picture on which our model is based on the splitting (or smearing) of the 
core line observed in the ordered phase would persist into the conducting phase, since 
the structural deformations of the lattice persist — the loss of order is a loss of three- 
dimensional coherence only and is essentially irrelevant to the site-local physics of core 
holes. 

Importantly however, a splitting (smearing) of the line due to the inequivalence 
of Ta sites will not alter the asymmetry of the underlying lineshape. // there are 
zero energy charge excitations in the conducting phase the measured split (smeared) 
lineshape should also be asymmetric. Shifts in the line due to changes in the local 
susceptibility will persist for the same reason. In bulk Ta the 4/ orbitals may be fitted 
with a Doniach-Sunjic lineshape and has an asymmetry aE ~ 0.1. 
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There are a number of other effects which need to be borne in mind. Foremost 
among these is spin-orbit couphng which can spht core hnes by large energies. In the 
case of the Ta^f state discussed above spin-orbit couphng leads to a splitting of ~ 3eV 
between 4/7/2 and 4/5/2 states. 

We should also be aware that an attractive potential (such as a core hole) in a one 
dimensional electron system will generally have a bound state. The presence of this 
bound state will modify the threshold for photoemission out of the core level (giving in 
general two thresholds, one for sites with the bound state filled and one for those where 
it is empty) and, in a truly metallic system, will also modify the asymmetry measured 
in each of these cases As splitting due to bound states will be uncorrelated with 
the lattice distortion, it can in principle be distinguished from one due to inequivalent 
sites in a distorted lattice. The possibility of a line split by a bound state appearing as 
a single (asymmetric) line in XPS carried out at finite resolution poses a more serious 
threat to the interpretation of data on these systems, and would need to be considered 
before very strong conclusions were drawn about any measured asymmetry. 



8. NMR 



NMR provides a local probe of spin excitations which is exactly complementary to XPS 
(a probe of local charge excitations). In fact, up to a structure factor F{q), the nuclear 
relaxation rate 1/Ti is given by precisely the same formula as we derived in [Appendix D 



for the XPS asymmetry, but with the spin susceptibihty Xa substituted for the charge 
susceptibility Xp the structure factor F{q) substituted for the core-hole conduction 



electron interaction matrix element | Vg 

„2 



1 _ r ,.^^^^..^{xAq,^)} 



, limV F(g) ^^"^^' ^^ (34) 

Here /^b is the Bohr magneton and the nuclear spin has gyromagnetic ratio 7. We shall 
assume purely on-site coupling. This leads to a momentum- independent coupling F. 
Then the rate is simply proportional to the slope of the absorptive part of the local spin 
susceptibility at very low frequencies. We shall neglect spin-orbit coupling so that the 
spin susceptibility is well-defined and isotropic. 

Since the single-particle Green's function for the LRA liquid is diagonal in spin 
indices, we may evaluate the spin susceptibility in terms of the density of states in the 
same way as we evaluate the charge suscpetibility in [Appendix D| to find 

jr^^ = 2nj'kBFnl. (35) 

where uq is the density of states at the Fermi energy. For a regular metal (Fermi Liquid) 
where no is not a function of temperature this gives a constant relaxation rate (Korringa 
Law). 
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The quantity Uq is however very strongly temperature dependent in the LRA theory 
(Figure |]) and so the temperature dependences of the relaxation rate and can be used 
directly to measure of the extent to which the pseudogap depresses the density of states 
at the Fermi energy. NMR measurements are in principle possible for (TaSe4)2l since 
5*6 at least is NMR active, but no data at present exist. 

The NMR Knight shift is determined not by the local, but by the uniform 
susceptibility 

AH 



rr X. (36) 

-no 

and so follows this in temperature dependence. This is plotted in Figure ^ for our 
paramatrization of the LRA model. 

In a LL, the temperature dependance of 1/TiT is dominated at low temperatures 
by the scahng dimension of the q 2kf part of the spin susceptibility, which is set by 
the anomalous dimension of the charge field Kp 

' lU^r'. (37) 



TiT V T , 

where Tq is the characteristic temperature at which the g ~ and q ^ 2kf susceptibilities 
are of equal size. 

Data for the quasi one-dimensional organic spin density wave (SDW) system 
(TMTSF)2CI04 have been interprited as showing LL behaviour over a range of 
temperatures just above the SDW transition temperature, with the rather small value 
of Kp ~ 0.15 found from a fit to 1/TiT. As far as we are aware no such analysis of 
NMR data has been attempted for a CDW system. 



9. Comparison with other Theories 

We have presented predictions for XPS and NMR. This has a particular goal. In a 
theory such as the LRA model the pseudogap is felt equally in both spin and charge 
channels, and should be evident in the suppression of both XPS asymmetry and nuclear 
spin relaxation. Importantly, both are local probes of the charge (spin) susceptibility 
and as such, directly comparable. 

In the limit of large coherence length the LRA spectral function (7) reduces to the 
sum of two lorentzians. In the language of many-body Green's functions we may say 
that the poles associated with fermionic quasi-particle excitations have migrated from 
the real axis into the complex plane, becoming purely imaginary for zero frequency. 
The fact there are still simple poles in the electron Green's function means that these 
quasiparticle excitations are electron-like in nature and therefore carry both spin and 
charge. If the effects of the mean field gap persist in the spin channel for the metallic 
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phase, they will also persist in the charge channel, and the (quasi-)gaps measured in each 
channel should have the same temperature dependence — in the language associated 
with high-temperature superconductivity, the "quasi-gap" and the "spin-gap" should 
be in exact correspondence. 

We can make this statement more formal by considering the ratios of the XPS 
assymetry and the NMR relaxation time 

a(T) , , 

R = — ^ 38 

Within the LRA model (or any other model in which the pseudogap is felt equally in 
spin and charge channels) R should be essentially temperature-independent. 

Almost none of these features are shared with the canonical one-dimensional model 
of interacting electrons which we may most meaningfully compare with experiments 
on (TaSe4)2l (and other quasi one- dimensional charge density wavesystems) — the 
Luther-Emery (LE) model. Here there are no fermionic quasiparticles and the gap 
exists uniquely in the spin excitations of the system. This means that the ratio R of 
Equation (38) should be very strongly temperature dependent, due to the activated 
behavior of spin excitations. 

A plausible form for the spectral function of the LE liquid capturing the essential 
features of the gap in spin excitations the has been proposed by Voit ^Ijl . The spectral 
function has two dispersing maxima, one associated with charge and the other with 
spin excitations, but now only the spin peak is a true singularity; the divergence at the 
charge peak is cut off by the spin gap. 

The spectral properties of a related field-theoretical model of electrons interacting 
with phonons in one-dimension have also been derived by Wiegmann [H2|. This goes 



somewhat beyond the considerations of the other theories alluded to in this article in 
identifying the readjustment of the order parameter which accompanies the introduction 
or removal of an individual electon. It is characterized by extremely broad and 
asymmetric dispersing features, even in the absence of electron-electron interaction. 
The extreme breadth of the dispersing features in this model resembles photoemission 
spectra taken on the high-Tc cuprates rather more than those taken on quasi one- 
dimensional conductors. 

Because of the broadness and asymmetry of features in the LE spectral functions 
and its dual peak structure, it would be difficult to extract a single energy scale from it 
empirically. An experimental spectrum taken at finite resolution might well "measure" 
a different gap scale and temperature dependence different from that determined by 
transport, optical conductivity, or susceptibility experiments. 

Thus the LRA model offers a better description of existing photoemission data for 
this system [Q. The LE model may prove to be more nearly applicable to the other 
most commonly studied inorganic quasi one-dimensional charge density wave system. 
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Blue Bronze (KMnO^) |]T2[, but faces some challenge in explaining the suppression of 



charge excitations over the energy range 0.1 0.3eV seen in optical conductivity. 
10. The gap discrepancy 

There remains the discrepancy in gap values mentioned at the outset. 

First note that only the photoemission value is truly out of line. Based on this we 
may put forward three speculative resolutions. 

One is that photoemission, which takes place on a very fast time scale and is a high 
energy process which may couple to many excitations besides the formation of a hole, 
cannot be compared directly with a long-time scale, low energy experiment such as 
transport. This point of view is broadly compatible with the treatment of the electro n- 
phonon problem performed by Brazovskii and Dzyaloshinskii |^3| in which there is a 



clear separation of time scales between different electron and phonon excitations of the 
system. 

Second, a movement of intensity away from the Fermi energy in photoemission 
spectra could originate in "extrinsic" energy loss processes, i. e. , those affecting the 
outgoing electron. In the range considered, Drude and phonon losses are the most 
important loss mechanisms. These extrinsic processes may be independently measured 
by electron loss spectroscopy . 



Third, midgap states may explain the mismatch between photoemission and 
transport gaps. These are highly localized states in the middle of the gap, which 
cannot carry a current and (due to their small number) have a small cross section 
for photoemission, but which offer a reservoir of charge that can be thermally excited 
over (half) the gap and participate in transport. This conjecture leads directly to the 
prediction that the activation energy measured at low temperatures is of order A and not 
2 A. The transport activation energy is indeed approximately half the gap measured in 
photoemission, and a possible mechanism for the generation of localized mid-gap states 
does exist in (TaSe4)2l through interaction with the lattice distortion. There are four 
equivalent ways of accomplishing the nearly commensurate tetramerization of the lattice 
which is observed experimentally, and it is natural to suppose that this "frozen" phonon 
is divided into domains with different phase by kinks (solitons) in the order parameter. 
That kinks in an order parameter can bind (fractional) charge is well known and has 
been extensively studied in the case of polyacetylene |Q , which has a half filled band. A 
numerical study of the formation of mid-gap states in the nearly commensurate quarter- 
filled case was performed by Machida and Nakano EH] , who find fractal structure in the 



density of states reminiscent of the Hofstader problem |46|; in fact both the motion of 
an electron in an incommensurate potential in one dimension and in a magnetic field at 
incommensurate Landau level can be described by the same almost Mathieu equation. 
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As we have been arguing that the lattice distortion persists in the normal state of 
(TaSe4)2l (a weak superlattice corresponding to the condensed TA Phonon mode is still 
visible in X-Ray experiments at 300K) we can further speculate that the midgap states 
will persist to these temperatures, leading to the observation of the smaller "gap" in 
room temperature optical conductivity. However, this hypothesis has the weakness that 
the states involved must have sufficient weight to dominate transport properties but not 
enough to show up in photoemission. 

Note that tunneling experiments on (TaSe4)2l would be extremely informative. This 
experiment also measures the density of states. Given that photoemission is in conflict 
with other data, a check is necessary. 

Thus the different gap sizes measured by different probes in quasi-one dimensional 
systems are hard to reconcile with the LRA model unless some additional mechanism is 
postulated a posteri to explain the discrepancy. They may be somewhat less surprising 
within a strong interaction scenario, but the different determination of (pseudo-)gap 
sizes and temperature dependence by different probes is not a problem of the exotic 
non-Fermi liquid conducting phase only, but one which needs to be addressed also in 
the context of the ordered phase. In the metallic phase however, XPS provides, at least 
in principle, a novel means of exploring whether (quasi-)gaps exist for both spin and 
charge excitations or only for spin, and might be particularly informative if taken in 
conjunction with NMR experiments. 



11. Comparison with Underdoped High-temperature Superconductors 

What light does (TaSe4)2l shed on the pseudogap in high-Tc superconductors ? 
Photoemission results on the two systems are the obvious place to start. In underdoped 
Bi2Sr2CaCu208+x above Tc the peak of the spectral function A{k,u) never reaches 
the Fermi energy, stopping about ~ lOmeV below the Fermi energy when k is along 



the (0, tt) direction Ths precise size of this pseudogap depends on doping 

and temperature and vanishes at a temperature T* > T^. This observation is 
qualitatively similar to that in (TaSe4)2l. There is, however, more than an order of 
magnitude difference in the size of the pseudogaps. Furthermore, the pseudogap in 
Bi2Sr2CaCu208+x is anisotropic with little or no pseudogap observed in the (vr, vr) 
direction. This possibility does not arise in (TaSe4)2l. The similarities go beyond this. 
The uniform magnetic susceptibility is suppressed in underdoped high-Tc systems 



and in (TaSe4)2l 0, showing that the pseudogap is present in the spin channel as 
well. The temperature dependence of the resistivities p{T) is superficially very different. 
In the underdoped high-Tc systems, p rises rapidly and monotonically between Tc 
and T*, though there are strong deviations from the linear behavior of p{T) which 
is characteristic of the optimally doped material. In (TaSe4)2l, p is roughly independent 
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of temperature in experiments, as we have seen. However, this may only reflect the 
fact that the accessible regime of temperatures does not reach up to Tmf- We would 
expect p to rise when T ~ Tmf, as seen in Fig. |[ Finally, we may compare optical 
conductivities a{uj). The most interesting high-Tc data from our point of view are 
actually c-axis conductivities which have no Drude peak at low frequencies These 
show a very clear pseudogap opening up as a function of temperature. In the a-b plane 
conductivity there is a Drude peak. This interacts in a subtle way with the suppression 
of the density of states and the evolution of the electron lifetime due to the pseudogap, 
rendering the interpretation of the data somewhat complicated. In (TaSe4)2l we are 
more fortunate in that the Drude width is considerably less than the pseudogap energy 
and the structures in cr{uj) are therefore well separated. Again, given the difference 
in parameters, it appears that the two systems are rather comparable. NMR and 
tunneling ||52| also show evidence for a pseudogap in underdoped high-Tc systems. These 
experiments have not been performed in (TaSe4)2l, but would be illuminating, as we 
have already stressed. 

The very strong similarity between the two systems suggests a common origin for the 
pseudogap behavior In (TaSe4)2l, it is clear that the properties of the pseudogap 
phase are essentially determined by the fluctuations remaining from the low-temperature 
phase. There is a very large difference between the actual critical temperature and 
the mean-field temperature. This is expected in this quasi-one-dimensional material. 
In the high-Tc systems, the difference is not so large, as the system is quasi-two- 
dimensional. If we take into account this dissimilarity it appears that in the intermediate 
pseudogap region, the fluctuations from the ordered phase are the determining factor. 
Thus, the comparison supports the point of view that the origin of the pseudogap in is 
superconducting fluctuations. The anisotropy in the pseudogap is consistent with this 
as well. 

Spin-charge separated theories of the pseudogap regime of high-Tc superconductors 
analgous to the LE liquid have been advanced by several authors fS^. The comparison 
of XPS and NMR data for these systems along the lines which we suggested for a one- 
dimensional charge density wave system might also shed light on whether the (pseudo- 
)gap exists in both spin and charge channels. We hope to develop this idea further in a 
later paper. 



12. Conclusions 

The Lee, Rice and Anderson model offers the simplest possible scenario for fluctuations 
of charge density wave order in quasi one- dimensional charge density wave systems that 
reduces simply to the mean field ordered state. The naive application of this model 
to existing data for the quasi one-dimensional charge density wave system (TaSe4)2l is 
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quite successful. In those cases where the LRA model fails to explain experimental data 
for this system above the charge density wave transition temperature {T^d), the usual 
mean field picture of the ordered state must also be called into question below T^d- 
The most important of these discrepancies is the size of the charge density wave gap 
(pseudogap) which is found to be of different size by different experimental techniques, 
both above and below T^d- 

This simple view of the origin of a pseudogap in a quasi-one-dimensional system 
should be contrasted with the canonical picture of "gapping" in 1-D, strongly correlated 
systems with strong CDW fluctuations, in which non-current conserving terms in the 
Hamiltonian lead to the formation of a gap for spin, but not for charge excitations. 
In this scenario, dispersing features in the electronic spectral function are generally 
extremely broad and asymmetric, and while the ordered state of the system must be 
three-dimensional, it need not be mean field like in nature. In those cases where we have 
been able to compare the LRA model directly to experiments on (TaSe4)2l it seems to 
offer a better description of data than existing predictions for strongly correlated models. 

XPS provides, at least in principle, a novel means of exploring which of these 
scenarios is correct, and might be particularly informative if taken in conjunction with 
NMR experiments. Tunneling measurements of the density of states would serve as a 
useful check on photoemission. 
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Figure captions 



Figure 1. a) Structuc of Ta chains in (TaSe4)2l. b) Schematic displacement of Ta 
atoms in tetramerized ground state. 



Figure 2. Example of ARPES spectrum for dispersing feature with k=ki? and 
theoretical fit, taken from Ref. 4. 



Figure 3. The temperature dependence of the uniform susceptibility as a function of 
temperature. The points are data and the line is the LRA model fit, both from Ref. 5. 



Figure 4. The DC resistivity over the experimentally accessible temperature range. 
The data in (a) are taken from Ref. 12 and the theoretical curve in (b) is plotted using 
Eq. 8. 
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Figure 5. a) Experimental spectra for the optical conductivity (TaSe4)2l taken from 
Ref. 8. b) The intrinsic optical conductivity of an LRA liquid as calculated from 
current-current correlations according to Eq. 12. 



Figure 6. The temperature dependence of the density of states at the Fermi energy 
squared normalized to free electron values, for the LRA model, as parametrized to 
describe optical conductivity measurements. The temperature dependence of NMR 
1/Ti and the XPS asymmetry ue, should both follow up to the accuracy of this 
parametrization of the model. 



Figure 7. Numerically determined change in the local static susceptibility 
J2q x(9) 0) T) relative to Tsd. This determines the shift in XPS lines. 



Figure 8. a) MP theory of the charge density wave state in terms of multiple 
scatterings between right and left moving electron states. These may be resummed 
to give a self energy correction which coincidentally is equivalent to b) the first term 
in a perturbation theory developed using the a; — > limit of the phonon propagator 
D{q) ~ Sq±2kfS{u}). 



Figure 9. Diagrams which must be evaluated when using the oj ^ limit of the 
phonon propagator, as in the theories due to Sadovskii. Propagators may be modified 
to take account of the finite correlation length of the lattice distortion, but the theory 
is then no longer exactly soluble. Crossing diagrams at higher order are included. 



Figure 10. Diagrams calculated in the "static non-crossing " approximation with the 
phonon propagator is Sq±2kfS{u}). 
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Appendix A. Formal standing of the LRA theory 

In this appendix we discuss the nature and vahdity of the LRA theory from a technical 
rather than phenomenological point of view. 

We begin by rederiving the self energy correction (6) using equations of motion in 
the spirit of the original LRA paper, but showing that on a more careful examination 
of the nature of the averaging involved, the equations of motion do in fact close. 

Our starting point is the Hamiltonian 

i^LRA= E + ^QC-fe'+QC-fe'] (Al) 

k Q,k'>0 

= -^g{Q)HQ))- (A2) 

and 

The imaginary time dependence of the single electron Green's operator (in the 
Heisenberg picture) is given by 

G{k\k,r)^ -T,Ckir)cl{0) (A3) 

= - {e(r)c,,(r)4(0) - e(-r)4(0)c,,(r)} (A4) 

and its evolution follows 

= - Hr){cA0)Am - r.^4(0). (A5) 

The Fermi field in turn has dynamics governed by 

Ck'ir) =e^^Cfe,e-^^ (A6) 
^^ = e^^[//,c,,]e-^^ (A7) 

and the phonon field, by exphcit assumption, is static — i. e. has no dynamics. 
Inserting the Hamiltonian (Al) in (A6) we find 

dGik', k, r) ^ _ ^^^^^^^^ _ (^g^ 



dr 



Y.^*QG(k' + Q,k,r) + h.c.\. (A9) 

Q 



Then for the Green's operator with k' — k. 
dG{k, r) 



dr 



= -S{t) -e{k)G{k,T) (AlO) 
Y,^*QG{k + Q,k,T) + h.c.\ (All) 
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We introduce the Fourier transform 



G{T) = -Y.G{iujr,)e-^-- (A12) 



Whence 



-iUnG{k,iuJn) = — 1 — e{k)G{k,iuJn) (A13) 

J2'^*QG{k + Q,k,iujr,) + h.c.\ (A14) 
Q ) 

Physically we are interested in scattering between the two Fermi points, i. e. k' ~ kf, 
Q ^ 2kf, i. e. 

[e{k) - iiUn]Gik, lujn) = 1 - E ^^^(A; - Q, k, iujn) (A15) 

Q 

Similarly, 

[e{k -Q)- iuJn]G{k - Q, k, iuj„) = - ^ %'G{k + Q' ~Q, k, iuj„) (A16) 

Q' 

which implies 

e(fc) - tuJn\G{k, tujn) = 1 + — n 7^^ ■ (^1^) 

e[k - Q) - tuJn 

We now turn this relationship between operators into an equation for the single 
electron Green's function by taking the expectation value of both sides 

[t[k) - lU}n\{{G{k,lUn))) = 1 + 2^ 7T Jy, ■ 1^8) 

QQi ^{i^ ~ Vj ''^n 

where ((...)) denotes both thermal and quantum mechanical averageing. 

In the case of the classical field \E'q quantum averageing is irrelevant and we are left 
with the thermal average (\I/q;^q)t over all possible configurations of the static phonon 
field. By definition {{G{k,iuJn))) = G{k,iujn) the temperature Green's function for the 
system. 

We note that for a many-chain system with "frozen" phonon disorder one might 
equally motivate an ensemble average on potentials from the fact that all experiments 
measure average over many chains, and the potentials on different chains would in 
general be different. 

Following Scalapino, Sears and Ferrell (SSF) |]55|, LRA consider consider a one- 
dimensional classical order parameter "^q described by a free energy 

FI^Hq] = a(T)|vl/Q|2 + 6(T)|v^q|^ + c(T)(Q - 2M'I^qI', (A19) 

where the parameters a(T), b{T) and c(T) are taken from the mean field treatment of 
the electron-phonon interaction. Using this free energy to perform the thermodynamic 
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average it is then found that correlations of lattice at order decay exponentially 
with distance 

(^(x)^(x'))r = {ip'^{T))cos[2kf{x - x')] exp[-|x - x'\C\T)] (A20) 

where the temperature dependence of the two parameters follows from {ip'^{T)) and 
^~^(T) may be found from an exact mapping onto a quantum mechanical anharmonic 
oscillator problem p5| . 

It follows from (A20) that 

and the equations of motion now close. We assume that within this thermal ensemble 
(\I/q) = which implies that there is no anomalous propagator above the three- 
dimensional ordering temperature T3D. 

The remaining sum over Q is performed as a contour integral, using the result that, 
for linearized dispersion 

t{k - 2kf) = - t{k) (A22) 

yielding 

lUJn + t{k) ±lVft, [I) 

with the choice of sign depending on the half plane to which the function is to be 
continued. 

The self-energy correction which we derived above 

^LRAik, tU^) = . ^ (^24) 

has the correct limits for ^^^(T) (BCS mean field theory with gap = {ip'^iT))) 
and for {ip^{T)) — > (free electron gas). Physically it can be thought of as describing 
the pairing of a coherent electron of momentum k ^ kj with an incoherent hole of 
momentum k ^ —kj . However since it at first sight a somewhat crude theory, our use 
of it requires some clarifcation and justification. 

The physical motivation is clear - below the mean field transition temperature 
there is a very strong separation of time scales between the soft phonon modes and 
the electrons. This observation might equally inspire us to develop a perturbation 
theory for electron-phonon coupling in which the electrons scattered elastically from 
the fluctuations of lattice order. This would require using the full electron-phonon 
vertex of the Frohlich Hamiltonian (1) and the u ^ limit of the phonon propagator. 
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Indeed if we interpret the result of the thermal average over static phonon fields 
A21 as an — phonon propagator 

D{Q) = mT)) ^~'^^\ o,\2^,-2 (A25) 

then the lowest order self energy correction for electrons scattering off this (quantum 
mechanical) phonon field is indeed the LRA result (A24), and the LRA theory is often 
discussed in these terms. 

It is not hard to calculate the corrections to this self-energy which occur at the next 
(and higher) order in this interaction They are not small, and tend to "wash out" 



the pseudogap found in the LRA theory — why then should we trust it ? 

Under the model assumptions made above about the thermal ensemble of lattice 
potentials, the equations of motion for the Hamiltonian (4) close. This means that the 
LRA theory we present is the correct a solution of a slightly different problem. We 
explain why, and what this problem is, below. 

By using the expectation value of the phonon field in the Hamiltonian (4) we are 
going somewhat further than taking the u —>■ limit of a quantum mechanical problem; 
we are removing the phonons from the quantum mechanics of the problem entirely. The 
"frozen" lattice is treated strictly as a static external field - a classical object. 

This approach has the great advantage of reducing exactly to the familiar "mean 
field" case in the fully ordered limit where 

= A5Q±2fc^. (A26) 

in which case the Hamiltonian is 

H=Y: <k)clcu + E[A*cL2fc,c, + Ac^+2.,c-fe] (A27) 

k k 

A = -l=g{2k^){u2kF). (A28) 
y L 

This familiar Hamiltonian may be solved exactly by the equation of motion technique 
used above, or by canonical transformation, and one obtains a self-energy of the form 



A 



|2 



Y.MF{k,iuJn) = -. ■ TTT (A29) 

lUn + e{k) 

The same self-energy correstion is indeed the result found in a perturbation theory 
for electrons interacting with a limu; — > phonon propagator 

D{Q) = A5Q±2fc/(^) (A30) 

but only at lowest order. Resumation of the remaining diagrams at higher order 
does not lead to recovery of the correct BCS-like mean field form for the electron Green's 
function. Instead one finds a propagator for a system with a distribution of real gaps. 

G{k,iun)= rdxe-^-^^ — (A31) 
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a result that first published by Sadovskii [^. A later attempt by Sadovskii ||58[ to 
generalize this treatment to the case with a finite coherence length (phonon propagator 
A25) has been widely discussed but seems to be fiawed |]59 |. 

The reason for the apparent contradiction between the exact result and this 
resumation of perturbation theory is that so far as Hamiltonian (A27) is concerned, 
the wrong perturbation theory has been resummed. The topology of the diagrams for 
scattering off an external field (in our case "frozen" phonons) and for scattering from 
an — * phonon propagator are different because a propagator (two points joined 
by a line) explicitly correlates pairs of scattering events, whilst all scatterings from an 
external field are independant of one another. 

In solving the mean field Hamiltonian we have, in diagrammatic langage, resummed 
the infinite series of all possible mulitiple scaterings which transfer electrons between the 
two Fermi points (see Equation (H)). The vertex in (A27) is of the form of an external 
field which scatters particles from the right Fermi point to the left (or vice versa) whilst 
preserving the new momentum q. 

If one wishes do perturbation theory starting from the Hamiltonian one should 
consider diagrams of the same type, for the same reason. The coincidence between the 
first order self energy correction a perturbation theory developed using an u; ^ phonon 
propagator with Lorentzian form and the LRA self energy is in some sense just that 
- a coincidence - and moreover the same coincidence as the correspondence between 
the exact self energy of the usual mean field picture and the lowest order self-energy 
correction for interaction with a phonon propagator of delta function form. 

Our goal is to understand the effects of fluctations of charge density wave order in 
a quasi-one dimensional system. In the absence of any exact solution to which we can 
compare, whether it is a better approximation scheme to remove quantum mechanics 
from the phonon field altogether and so be able to recover mean field theory, or to 
consider a problem in which both electrons and phonons are quantum mechanical objects 
but scattering between them is elastic remains an empirical question. As the thermal 
averaging over static field configurations in the LRA scheme in some ways imitates the 
averaging over many chains which must take place in the real system, even an exact 
solution which called the LRA scheme into question in ID might not invalidate its 
application to real systems. 

A comprehensive treatment of the Sadovskii model and related technical issues is 
supplied by Tchernyshyov |^|. A slightly different approach to the same problem has 
also been considered by Millis and Monien ||60|| . 
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Appendix B. Parametrization of Model 

The parametrization of the model involves determining the temperature dependance of 
the two parameters of the self energy correction, {tp'^{T)) and the coherence length 
^~^(T). This is to some extent arbritrary. As we intend to use the model in a 
phenomenological but self-consistent way they should ideally be taken from experiment; 
in order to have a simple working picture we borrow the analysis presented in the original 
LRA paper which leans heavily on the results of SSF |^ . 



The classical field \1/q has a free energy with parameters taken from the mean field 
(linear response) perturbative treatment of the ID Frohlich Hamiltonian : 



with 



F[^q] = a{T)\^Q\' + + c{T){Q - 2^1^^!', (B32) 



T — T 
a{T) = D,-—^ 



KT)=Do[bo + ibo-h)^ 



c 



and 



c{T) = DoeoiT), 



where Dq is the (constant) density of states for the band, which is taken to have width 
2eF- The parameter ^o(^) is the length scale emerging naturally from the linear response 
analysis of the one-dimensional electron phonon problem 

eo(T) = ^ (B33) 



We fix bo and bi to give the correct zero temperature value of the gap Aq for the 
experiment under consideration (A^-'^^-^'^ 7^ ^transport^ according to 

60 = ^, (B34) 

and 

7C(3) (1.76)2 



61 = bn- 

1671 0.5 

The problem of determining {ip'^{T)) and ^^^(T) then reduces to that of finding the 
low lying energy levels of a particle moving in an anharmonic potential well, the shape 
of which is determined by the coefficients of the free energy 

H = - + a{Tm' + b{Tm' (B35) 



34 



This problem may solved numerically, or approximately using perturbation theory and 
assymptotic analysis. 

We consider following parametrization, found from a simple perturbation theory 
sufficently accurate in the temperature range of experimental interest in the next section: 

4T 1 

r'(r)=f„-'(r)(— --) (B36) 

where a' = a{T)/T. These approximate forms of the parameters {tp'^{T)) and ^^^(T) 
are used wherever we compare with expermiment. 

One interesting feature of the exact (numerical) determination of the ^(T) is that 
the correlation length is strictly infinite only at T = but begins to diverge strongly 
at a temperature approximately one quarter of the mean field transition temperature. 
This fact leads H to the simple prediction that 



- 3D 



Tmf/4. (B38) 



The value of 263 K for Tgo found by experiment (Table 1) is therefore compatible within 
a fluctuating charge density wave scenario with the estimated Tmf of ~ 900K. 

An interesting associated technical question is whether, in a rigorous treatment of 
the model defined by Equation (1), ^(T) actually does diverge at zero temperature, 
or whether it remains finite due to quantum fluctuations. This depends on properly 
including the Umklapp processes present in a quarter-filled band. However, for the 
temperature range T > T^d which is of interest in this paper quantum fluctuations are 
not important. For T < T^u, three-dimensional effects dominate in the real system, 
again suppressing quantum fluctuations. 

Appendix C. Static noncrossing approximation 

Another possible way of calculating the single-particle Green's function in the static 
limit is to replace the phonon propagator with delta functions: 

D{q,u;) = 5g±2kJiu;). (C39) 

and in addition to make the noncrossing approximation that no phonon lines cross. 
This amounts to calculating the Green's function self-consistently but neglecting vertex 



corrections entirely. See Figure 10 for an illustration of this. We shall find that this 
method leads to unphysical results. It is mentioned only because it is one way of 
obtaining a pseudogap in a number of different models. 
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The self energies for two wavectors are 

UK^) = — ^1 — 7 (C40) 

- ek-2kF - ^{k - 2A;f, uj) 
^(k - 2kF, UJ) = ^-T^Ti r. (C41) 

Because of the delta-function effective interaction, these equations close. 
Consider first the special case k — kp- Then I^{k,u!) — E(A; — 2kF,uj) — E(— /cp, a;), 
ek — and 

^kF,u;) = ^ (C42) 

UJ — h[kF, UJ) 

which is a quadratic equation with solution 

j:{kF,u;) ^V^2_4^2. (C43) 

The self-energy has a finite imaginary part for |a;| < 2|gi|, i. e. , even at zero frequency. 

The general case is only slightly more involved. Let us define q — k — kp so that 
Cfe = vpq €k-2kF = —vpq- The result is: 

1/2 

2 / a; — Vpq 



1 1 



{u - Vpq) - 4:g 



UJ + Vpq 



(C44) 



The imaginary part of this function is nonzero when \uj + Vpq\ < 2g. 

A finite relaxation time makes no sense for a static, ordered, potential. We conclude 
that the noncrossing approximation is not appropriate for this problem. This is not 
surprising. Neglect of vertex corrections is only valid when the electron wavefunctions 
are only slightly perturbed by the external potential. When degeneracies {k — kp) or 
near-degeneracies {k ~ kp) are a dominant effect, ordinary nondegenerate perturbation 
theory is not accurate. 

Appendix D. XPS Lineshape and Shift 

We define a core retarded core hole Green's function 

Gh{t) = -ie{t){d\t)d) (D45) 

and divide it into two parts 

Gh{t)^ -ie{t)e-^^'p{t) (D46) 

where ujt is a threshold energy (containing, in general, contributions from all orders 
of perturbation theory), and all the relevant many particle effects are consigned to the 
factor 



p{t) = (I Texp 



i f dtiViti) 
Jo 



(D47) 
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The lowest order term in a linked cluster expansion for p{t) is 



Pit) 



Ht) = U-^f fdh fdh{TV{h)V{h)) 
Z Jo Jo 

= - ^ E I r (Tpiq, U - t2)p(-q, t2 - h)) 



(D48) 
(D49) 



We will confine ourselves here to examining the contribution to the Fermi Edge 
Singularity (FES) from this term, which contains all the essential physics of the problem. 

The expectation value {Tp{q,ti — t2)p{—q,t2 — ti)), obviously has the form of a 
density-density correlation function. We therefore define 

ill{q, iVLm) = - T r dre'^'-^iTpiq, t, - t2)p{-q, ^2 - ^i)) (D50) 
L Jo 

and introduce a bosonic spectral function B{q,uj) 



i-r/ \ dw 



Biq,U})_ ^ B{q, -u) 



(D51) 



where iflm is a bosonic Matsubara frequency. 

Using the general result B{q,u;) — —B{q, —Cu), ( which follows from analyticity and 
parity), we find that 



pt rt poo 

/ dti dt2U{q,ti-t2) ^2 dCjB{q, 
Jo Jo Jo 



-iuit 



CO 



and therefore 



1 roo 

F2{t)= -yElKl' / dCuB{q, 



1-e 



-iCjt 



(D52) 
(D53) 



where we have dropped terms linear in t (and therefore belonging in ujt). 

In order to be able to apply our analysis to a system with arbitrary spectral function 
A{k^ oo) we rewrite the correlation function n(g, iVLm) in terms of Fermionic propagators, 
and then substitute a spectral representation of the electrons 

duj' ' A{k,uj') _^ A{k, -oj'' " 
10 27r iujn — LO' iujn + oj' 



(D54) 



we find 



/CO (j^' 
—A{k + q,u;' + u;)A{k,u;') 
-oo ZTT 



X [nf{uj')-nf{u' + u)] 

We Taylor expand 

E b^fi^k+q) - nf{ek)] 6{uj - Sk+q + efc) 

kq 

= -^E -^^(^ - ^k+q + efe) + . . . 
kq 



(D55) 



(D56) 
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At T — 0, duf/duj — —5{ijj)] and we can eliminate both integrals 

B{q,uj) = ujY,A{k + q,uj)A{k,Q) (D57) 

k 

We now make the assumption that Vq has no significant q dependence and define 
the function 

BM=E^ (D58) 

q ^ 

In general B{ijj) will also be a function of Z, where I is an index over scattering channels 
(in ID the limiting values q ^ 0, q ^ 2kf). For simplicity we suppress this dependence 
in what follows. 

In terms of B we may write 

F2(t) - - I K) P / duB{u) z (D59) 

L Jo oj 

In the physically relevant limit (a; ^ 0,i ^ oo), and making the substitution 

k' — k + q we find 

1 2 



j:A{k,o) 

. k 



nl (D60) 



B{u; ^0)^J2 Mk', OJ 0)A(A;, 0) 

kk' 

where Uq is the exact density of states at the Fermi energy. 

We impose a soft bandwidth cutofi^ eo (equivalent to having made the substitution 
I K) P e~''/'^° K) P e~'^/'^° in D48), in which case the integral over a) can be 

performed exactly to give 

F2{t) = - I Vb I' ln(l + ieot) (D61) 

Therefore, 

G.(t)- -ze(t)e--*^-L_ (D62) 

where for a ID electron system within these approximations and restoring the 2k f 
scattering channel, the exponent is given by 

ocE^nl{\V^\'' + \V2k,?) (D63) 

In the special case of free electrons we see that we have recovered the Born 
Approximation to the exact result. 
More generally 

a,= limlV|y,r^M^ (D64) 

provided that a) the unperturbed Hamiltonian is quadratic in Fermi fields b) the leading 
small uj behaviour of Y,q I ^ P '^{xpio.-i^)} is linear. 
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This form of Gh{t) leads directly to the power law assymetries in core hole hneshapes 
observed in experiment. In the case of non-interacting electrons it can be shown that 
second order perturbation theory represents the Born approximation to the the result 
obtained from resuming all orders of perturbation theory. 



TT 



(D65) 



where Si{Ef) is the exact phase shift at the Fermi surface in the Vth scattering channel. 

It is also possible to calculate the second order shift in the XPS line due to its 
interaction with the conduction electrons. This may be evaluated as a conventional self 
energy correction; the result is 



(D66) 



Gross changes in the susceptibilty due (for example) to the opening of a gap in a 
previously ungapped system, can therefore lead to a shift in core lines as well as a 
change in their assymetry. 



Appendix E. Core Level coupled to a Luttinger Liquid 

The following model for a localized single Fermion d (chosen here to be a hole) interacting 
with a set of Bosons a„ 



H 



ee + $]Mg(a, + aL,) 



may be solved exactly by the canonical transformation 
H = e'He-' = dd} [e^ - A] + Y^^A^i 



where 



(E67) 

(E68) 

(E69) 
(E70) 



This model has the unphysical feature that the interaction term does not conserve 
momentum - the core level does not recoil when electrons are scattered from it. 
It is also possible to find the fermion correlation function 

iGdit) = {Td^it)d) (E71) 

in closed form. We consider only the case t > 0, for which the result is 

iGd{t) = [1 - n/(ee - A)] exp [-(f){t)] 



(E72) 
(E73) 
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where Nq is a boson occupation factor, which we set to zero in what follows. The 
procedure for solving the model and obtaining correlation functions is described in some 
detail in |61|. 

The simplest prototypical Hamiltonian for a core level coupled to interacting 
conduction electrons in ID may be written 



(E74) 

(E75) 
(E76) 



where ep is the free electron dispersion and Vg the FT of the electron-electron interaction. 
This model (the Tomonaga Model) may be solved by linearizing the electron dispersion 
ep = Vf{\p\ — pf) and introducing a set of collective Bosonic coordinates to describe 
excitations of the conduction electrons. 

We now make the observation that the bosonized form of this model is (up to a 
constant) in exact correspondence with model of a single fermion interacting with a set 
of bosons solved above. 

Explicitly, we split the density operator into separate parts referring to right and 
left moving electrons p{q) = pi{q) + P2(?) which we normalize as bosonic operators 



TT 



-r[O{q)pi{q) + 0{-q)p2{q)] 
Ltt 

6. 



qq 



We write the kinetic energy term as 



(E77) 
(E78) 

(E79) 



and make use of the canonical transformation 



to obtain 



H=Y.Eq 



lUJa 



' — {ag + a'_gj 

Eg 



ag, a'q, 



Vf I q 



^qq 



J2Mg{ag + alg) 



+ eddd'^ 



(E80) 
(E81) 

(E82) 



which is clearly of the same form as (E67). 

The density of states at the Fermi energy in a Luttinger Liquid is identically zero and 
one might expect, on the basis of the arguments presented in Appendix D that this would 
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mean that it would present no asymmetry in core hole responses. This expectation is 
not born out — infact there is no essential difference between the bosonized Hamiltonian 
for free electrons coupled to a core level and that for interacting electrons coupled to 
a core level (provided that we neglect all "backscattering" processes). The response of 
the interacting system is therefore qualitatively the same. However, since the Bosonic 
eigenstates of the interacting electron system are not exactly the same as those for the 
free one, we must remember to modify the energy levels and matrix elements in the 
expression for Gd{t) E72. 

(E83) 
(E84) 

The all important logarithm in the response of the electron system to the impurity 
may be obtained by replacing the remaining sum over q with an integral. In the case 
where all interactions are S functions it is particularly easy to understand what happens. 

Then, 

Eq =v\q\ Mq = ^ (E86) 

Zn 

where v = Vf^l + 2V/'n'Vf and go is a band cutoff. We may bring this into the from of the 
non-interacting case by the change of variables q = {vf/v)q. Then, without needing to 
evaluate the integral, we know that the asymmetry exponent will be modified according 
to 

MV) = f^V«i^(0) (EST) 

3 

(E88) 

V / \ TTVf ) 

i.e. suppression for V > and enhancement for < 0. 

A more careful analysis of the g — > limit of integrand in the general case leads to 
the result (31) in terms of the Luttinger Liquid correlation parameter of the quoted in 
the text. The problem of the X-ray response of a LL has been considered previously by 
a number of authors ^ . 



V 



It is clear that the approach of Appendix 1 fails dramatically and it is worth pausing 
to consider why. The result (D53) is quite generally valid at this order. The error 
arises in the evaluating density-density correlation function n(g, iQ) as a single loop of 
fermions. In the case of a system with electron-like quasiparticles this is unlikely to lead 
to a large qualitative error, but in the case of a Luttinger liquid where Wick's Theorem 
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does not hold it is quite simply the wrong thing to do. Fortunately it is easy to evaluate 
the exact form of Il{q, i^) by bosonization. Substituting this in (D53) leads to the same 
results as the approach outlined above. 

We note that none of these arguments are affected by the introduction of spin, and 
may even be applied to an "insulating" phases, if the gap is in the spin rather than the 
charge sector (Luther-Emery liquid). 
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